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We generalize the cosmic energy equation to the case when massive particles in- 
teract via a modified gravitational potential of the form <p(a,r), which is allowed to 



explicitly depend on the cosmological time through the expansion factor a(t). Using 
the nonrelativistic approximation for particle dynamics, we derive the equation for 
the cosmological expansion which has the form of the Friedmann equation with a 
' renormalized gravitational constant. The generalized Layzer-Irvine cosmic energy 

r S^ , equation and the associated cosmic virial theorem are applied to some recently pro- 

^ ■ posed modifications of the Newtonian gravitational interaction between dark-matter 

£3 ' particles. We also draw attention to the possibility that the cosmic energy equation 

may be used to probe the expansion history of the universe thereby throwing light 



on the nature of dark matter and dark energy. 
PACS numbers: PACS number(s): 04.50.Kd, 98.80. Jk 
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I. INTRODUCTION 

Almost half a century ago, Layzer & Irvine [1], and independently Zeldovich & Dmitriev 
demonstrated that in an expanding universe, the peculiar kinetic (K) and potential (U) 
energies of a large system of pressureless particles interacting via the Newtonian potential 
= —G/r satisfy the cosmic energy equation 

E=-(2K + U)H, E = K + U, (1) 

where (see also 0,Q|) 

K=\Ys m ^' ( 2 ) 



G f [p(r 1 ) - g] [p(r 2 ) - g) , , 



and 

2 J |ri - r 2 \ 

are, respectively, the peculiar kinetic and potential energies of a system of particles rrii having 
coordinates and peculiar velocities Vi. Here, p(r) is the mass density of these particles, 
q is its cosmo logical background average, and H = a/a is the Hubble parameter describing 
the cosmic expansion. For the discrete distribution p(r) = Yli m i^( r ~ r i)i the integral in 
(EJ) should avoid the configuration T\ = r%. 

One can introduce the correlation function £(r): 

([p(r + r)-g][p(r )-g}) = g 2 ar) (4) 

with the obvious property 

J ^( r ) d 3 r = An J £(r) r 2 dr = . (5) 
Then Eq. ([T]) can be averaged to yield the corresponding equation per unit mass 

(£) = -[2{JC) + {U)]H (6) 

with 8 = E/M, U = U/M, and 

(U) = -2irGg J £{r)rdr . (7) 

Here, M = gV is the total mass of the system in a large volume V. 
Two important limiting cases of Eq. ([1]) need mention {3, ^]: 

(i) For noninteracting particles, we have {7 = 0, and the relation K = —2HK corresponds 
to the kinematic decay of the peculiar velocities with time Vi oc l/a(t). 
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(ii) Once the clustering has entered a stationary regime, the condition E = results in the 
virial relation 

2K+U = 0. (8) 

A breakdown of (JE]) for galaxies belonging to the Coma cluster led Zwicky to suggest 
that a large amount of dark matter might dominate the dynamics of Coma. While dark 
matter appears to be ubiquitous, its nature remains elusive. Indeed it is now believed that, 
to properly account for the energy budget of the universe, dark matter (DM) must be supple- 
mented by an even more enigmatic 'substance' called dark energy (DE), which, on account 
of its large negative pressure, causes the universe to accelerate instead of decelerating. 

Our current theoretical understanding of DM and DE can be broadly divided into ideas 
that are mainstream and those that are radical. Mainstream notions suggest that DM is 
formed of nonbaryonic particles of relic origin while DE consists of the cosmological constant 
or a relic scalar field such as quintessence. Radical notions suggest that the purported 
existence of DM/DE may be pointing to a breakdown of Newtonian/Einsteinian gravity on 
large scales. Theoretical models which incorporate this latter set of ideas include Braneworld 
and f(R) gravity theories in the case of DE, and Modified Newtonian Dynamics (MOND) 
and the screened gravitational interaction model in the case of DM (see also 

Ever since its discovery, the cosmic energy equation ([1]) has been one of the bulwarks of 
modern cosmology, and its many applications include estimates of the matter density and its 
gravitational binding energy 0,0,1s]. In this paper, we show how the cosmic energy equation 
can be generalized to incorporate more flexible forms of the gravitational interaction of dark 
matter, some of which have been suggested in the literature. 

Our assumption is that nonrelativistic dark-matter particles interact with each other via 
the two-particle potential 0(a, r) so that the potential energy between two particles is 

mim 2 (a, \r 1 - r 2 |) . (9) 

As indicated in (jHJ), we allow the two-particle potential to depend explicitly on the cosmo- 
logical time through the scale factor a(t), which is characteristic of some class of the models 
to be considered below. We consider the universe dominated by these dark-matter particles 
and by the uniform dark energy. 

In Sec. [Til using the nonrelativistic theory, we obtain the law of cosmological expan- 
sion for our universe, which turns out to have the Friedmannian form with a renormalized 
gravitational constant. In Sec. IHIl we then derive the generalized cosmic energy equation 
and the corresponding virial relations, which are analogs of the corresponding Layzer-Irvine 
equation (JT]) and virial relation (jSj). In Sec. HVl we comment on the results obtained. 
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II. COSMOLOGICAL EXPANSION 



In this section, we will derive the law of cosmological expansion given the interaction 
potential fl9]). The universe model can be taken to be spatially flat since we are dealing with 
spatial scales much smaller than the Hubble length. We describe particle dynamics within 
a large volume using the nonrelativistic approximation with respect to comoving particle 
velocities. 



A. Expansion law in the ACDM model 

The Lagrangian for a nonrelativistic particle is given by the expression {^J 

£ = l m f 2 -m$(r,t), (10) 

where $(r, t) is the full potential acting on the particle. In the ACDM (cosmological con- 
stant A + cold dark matter) model, in which dark-matter particles interact via Newtonian 
potential, for a pointlike matter distribution p(r,t) = ^2 i rrii5(r — 1%), the potential $(r,t) 
is given by 

•WH-^ti-f' (id 

3 3 

Proceeding to the comoving coordinates x = r/a and making a canonical transformation 

£ ->■ £ - 4- (-maax 2 ) , (12) 
one transforms the one-particle Lagrangian ( flOl) to 

C = -ma 2 x 2 — imp , (13) 

where the new potential is given by 

1 



y? = $ + ^a'ctx 2 . (14) 



In the Newtonian case, with account taken of pip , the new potential is 

G sr^ m, If.. Aa 2 



3 



V = —Y-^- + - aa-^)x 2 . (15) 

Y a ^ x-xA 2 V 3 / 



The law of cosmological expansion can be derived directly from ( 1151) under the require- 
ment that the potential tp should not exert force on a test particle in the limit of a homoge- 
neous distribution of matter p(r) = g. To see this, we note that, in particular, the Laplacian 



5 



of the potential should vanish in this limit. Now, we have 

VV = AnGa 2 rrij 5(r - r,-) + 3ad - Aa 2 

j 

-)■ AnGa 2 g + 3ad - Aa 2 . (16) 

Whence, we obtain one of the Friedmann laws of expansion 

a A-kG A 

- = — ^ + ( 17 ) 
a 6 6 

B. Expansion law in the case of general potential 

In the case of a general two-particle potential the derivation that led to the cosmo- 
logical potential f|T5|) for a particle can be repeated and generalized to 

ip = ^m j( f)(a,a\x-x j \) + ^ (aa-^^-\x 2 , (18) 
j ' 
where A e ff is the time- dependent ingredient in the universe which describes dark energy and 
replaces the cosmological constant and, in particular, whose interaction with dark matter 
makes the potential (Q time-dependent. Our equations below will be general and will not 
require the specific knowledge of this term; we will only assume that it does not cluster, 
remaining homogeneous in space. One such field-theoretic example is considered in Sec. Ill CI 
below (example [aj) . 

The class of theories under consideration in this paper could, in fact, be defined by the 
property that the potential acting on a single dark- matter particle has the form ( fl8l) . This 
complicated potential already includes the response of the universe, with all its field-theoretic 
ingredients, to the presence of dark-matter particles at specified spatial positions. In this 
subsection, we are going to derive the modified expansion law by taking the Laplacian of 
(fT8l) and demanding that it vanish in the case of a spatially uniform distribution of matter. 

Suppose that the gravitational two-particle potential has the form (we suppress the de- 
pendence on a as it is not relevant here) 

0( r ) = _^ /(r) (19) 
r 

with a sufficiently regular function f(r) [in particular, lim rf'(r) = 0], as will be the case in 

i — 

all our concrete examples. Then, calculating its Laplacian, we have 

V^(r) = -Gf{r)V 2 r - - 2GV r - ■ V r /(r) - -V?/(r) 



2G G 
AnGf(0)5(r) + —f'(r)-- 



f"(r) + -f'(r) 
r 



AnGf(0) 5(r) - jf"(r) . (20) 
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Applying this to (1TB]) and requiring that V 2 <p = for uniformly distributed matter, we 
obtain the equation 

a 4:7iG cS A eff 

a = - — Q + — > (21) 
where the effective gravitational constant is given by 

G cS = G lim [/(r) - rf'(r)} . (22) 

The derivation works as long as the limit in (12 2 p exists and is finite, which we assume 
to be the case. For the Newtonian potential <fi = —G/r, we have / = 1, and Eq. (122"]) 
returns the usual gravitational constant G. However, if the potential decays faster than 
1/r at large distances, Eq. (I22p may imply G e s = 0, so that the gravity of matter will not 
influence the rate of cosmic expansion. Note that the effective gravitational constant will be 
time- dependent if the function / in flT9]) depends on time. 

C. Examples 

To give examples to which our results can be applied, we consider some simple modifica- 
tions of the gravitational interaction between dark-matter particles: 



(a) The screened potential discussed in [9| has the form 

0(a,r) = -^(l + /3e-^) , (23) 

where is a dimensionless constant of order unity, and the time-dependent screening 
length r s (t) is a comoving constant, i.e., r s (t) oc a(t), such that r s (t Q ) ~ 1 Mpc. The 
potential arises in the field-theoretic model of interaction of dark matter with dark energy 
via the scalar field [lo| in the version with two dark-matter families. A subdominant 
relativistic family is used to stabilize the value of the scalar field; then the dominant 
nonrelativistic dark-matter particles have constant mass and interact via gravity as well 
as via the scalar field so that the two-particle interaction potential ( 123]) is generated. 
A dark- matter particle then moves in potential (118j) with (p given by ( 123]) . so that our 
analysis is applicable to this situation. One arrives at the Friedmann equation ([2T]) with 
constant dark-energy term A e g, just as in the field-theoretic approach of 0, llo| . The 
potential ( 123]) leads to faster evacuation of matter from voids and to an earlier epoch 
of structure formation, which could be perceived as an advantage for this model over 
ACDM 

Using the potential (118j) and the equations of motion for a dark-matter particle x(t), 

x + 2Hx = —zV<p, (24) 
a z 



and proceeding along the same lines as in [3], Sec. 27, one easily obtains the exact nonlin- 
ear evolution equation for the Fourier components 5k of the density contrast p(x) / g — 1 
in this theory: 



4 + 2H8 k = ItcGq 



1 + 







$k + Ak — C k , 



1 + (a/kr s y 

where the nonlinear term Ay. and the velocity term C k are given by 



A k = AnGg ^ 
fcVo.fc 



kk' 



+ 



f3kk' 



k' 2 k' 2 + (a/r s y 



and 



(25) 



(26) 



(27) 



respectively. The linear part of ( 12"5"|) reproduces equation (6) of [9( and leads to a more 
rapid development of gravitational instability than in ACDM. 



(b) A power-law correction to the Newtonian potential on large scales 

r 



G 

r 



r + r 



n > 1, 



(28) 



which we have regularized on small scales to avoid a singularity. Potentials of this 
form with n = 2 arise in the Randall-Sundrum model [11] with a single large extra 
dimension. (However, in this case, r C 1 mm, making the correction on cosmological 
scales extremely small.) 



(c) A logarithmic correction to the Newtonian potential 



G faG , 

— log — | e 



■r/r c 



(29) 



The influence of the 'regularizing' exponent is confined to very large scales r c ^> r$. 



For all these potentials the limit ( 1221) gives 



GeS — G , 



(30) 



implying that the background cosmological evolution remains unmodified. 
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III. GENERALIZED LAYZER IRVINE EQUATION 



The Lagrangian of the many particle system is obtained by the summation of (fl3l) with 
account taken of ( |T8l) . ( 12TT) and ( 122]) and is equal to L = K — U, where the peculiar kinetic 
energy is [3, 4] 



K 



With consideration of Eq. (I2ip . potential (lisp can be written in the form 



V= I [ P {r')- Q }<P{a,\r-r'\)d 3 r' 



using which, one obtains 



r = ^ / [P( r i) - Q] [p( r 2) ~ Q] <f>(a,ru) d 3 r 1 d 3 r 2 



(31) 



(32) 



(33) 



where, for the discrete distribution p(r) = X]j m «^( r ~ r i)i the integral should avoid the 
configuration r x = r 2 . Here and below, we use the notation ry = |rj— Vj \ and = \xi—Xj\. 
The Hamiltonian of our system as a function of the canonical variables (x^p^ is given 

by 

n = K + U = ^J2^.+lJ \PM ~ 6] \PM - Q] «V(«, ax 12 ) d 3 Xl d 3 x 2 . (34) 

Using the fact that the product a 3 \p{r) — g] — ^2 i rriiS(x — Xj) — a 3 g does not depend 
explicitly on time, we immediately obtain the equation for the peculiar energy of the system: 



-2HK + j J Wi)-Q\W*)-Q\ 



{a,r 12 ) 



rn + 



{a, r 12 ) 



d 3 r 1 d 3 r 2 . (35) 



dr " da 

Similarly to (JHJ) and <$B), we can write the averaged equations for quantities per unit mass: 



(8) = -2H(K) + 2ixHq / i{r 



(a,r) d(f)(a,r) 
r H a 



dr 



da 



r dr . 



(36) 



(U) 



2M 



£{\ri - r 2 \)(f){a, \n - r 2 \) d 3 ri d 3 r 2 = 2txq / £(r)0(a, r) r 2 rfr . (37) 



Equations (j3"5"]) - (l3"TI) form the main results of this letter. 

In the Newtonian case, we have d<p/dr = — <p/r, d<p/da = 0, and Eqs. ( 135]) and ( 136]) 
reduce to the corresponding Layzer-Irvine Eqs. (0Q) and (EJ). In the general case, Eqs. ( 133]) 
and ( 136]) contain an unknown function <p(a,r). If one has a theory predicting the shape of 
4>, one can, in principle, use (135]) and ( 136]) to test it. 
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Once the system has decoupled from the Hubble expansion, its peculiar energy evolves 
mainly because of the time-dependence of the potential 0(a,r): 

(£) « 2nHg [ tir) 9 *^ ar 2 dr . (38) 
J oa 

Taking into account (136|) . we obtain the generalized virial relation in the form 

{K)=v Q Jt{r) 9 -^f±rHr. (39) 

Next, we apply the generalized cosmic energy equation to the modifications of the gravi- 
tational interaction between dark-matter particles that we listed in Sec. Ill CI 

(a) Substituting the screened potential (1231 into (1351) . one finds, quite remarkably, that 
the generalized cosmic energy equation reduces to its Layzer-Irvine form ([I]) with the 
modified total potential U determined by f )23|) and f )33|) . For a system decoupled from 
the Hubble expansion, the energy evolution, according to (1381) . is given by 



(£) = Jl^A I ^ r )e'^r 2 dr , (40) 



r 

and the generalized virial relation ( 1391) in this case reads 

2{K) + (U) = f ^ r )e- r/r °r 2 dr , (41) 

On length scales r C r s (t) < 1 Mpc, we recover the usual Newtonian relations: 

-!&- = 2(lQ+ (W)Newton = / ^dr = , (42) 

H r s J 

where, in the last equality, we have used property (JSJ), and 

(W)Ncwton = -2tt(1 + P)Gq J i{r)rdr (43) 

is the Newtonian virial gravitational energy per unit mass with renormalized gravita- 
tional coupling. 

It should be noted that, since r c ~ 1 Mpc roughly corresponds to the Abell radius 
associated with a cluster of galaxies, it is unlikely that fj42|) can be applied to galaxy 
clusters. On these scales, the full form of the virial relation (14"TT) should be used. 



(b) The power-law correction to the Newtonian potential ( |28l) leads to the cosmic energy 
equation: 

(6) w - \2{K) + (W) Ncwton ] H + 2tt(1 + n)GH 6 J £{r) {y^j rdr ' ( 44 ) 
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and to the virial relation 



2{K) + (W)Newton = 2tt(1 + n)GHg J £(r) rdr , (45) 

where (U) Newton is the usual averaged Newtonian peculiar potential energy per unit mass 
given by (JTj). 

(c) The logarithmic correction to the Newtonian potential (129]) leads to the usual cosmic 
energy equation 

{£) « - [2(/C> + (W) Newton ] H + 27TaGH P f ^ r )r 2 e- r ^dr « - [2(/C) + (W) N cwton] # , 

?"0 J 

(46) 

where the last equality is valid in view of (JH]) if the correlation length of the system is 
much smaller than the cutoff scale r r . The virial relation in this case is 



2</C> + (U) Newton — 



j i{r)r 2 e- r/rc dr w . (47) 



r 

IV. DISCUSSION 

We have investigated a cosmological theory which produces a modified and, perhaps, 
time-dependent gravitational potential between matter particles in the form (Q. For such 
a theory, we have determined the background cosmological equation, which turned out to 
have the Friedmann form possibly with a modified gravitational constant. We have also 
derived the cosmic energy equation generalizing the Layzer-Irvine equation for the theory 
under investigation. 

We also note that the cosmic energy equation (jBJ) or its generalization ( |36l) can be used 
to determine the cosmic expansion history. For instance, from ([1]) we obtain 

The Hubble parameter H(z) determined in this manner could shed light on the nature 

m 

of dark energy either through the Om diagnostic [12j or by means of the effective equation 
of state of dark energy [l^ 



, x (2x/3)d\nH/dx - 1 

^''- i- ww ' x = 1 + z - ( ) 

This determination of the cosmic expansion history H(z), via the energy equation, is com- 
plementary to usual methods which rely either on standard candles (supernovae of type la) 
or rulers (baryon acoustic oscillations). This possibility will be examined in detail elsewhere. 



11 



ACKNOWLEDGMENTS 

The work of Y. S. was supported by the Cosmomicrophysics programme and by the 
Program of Fundamental Research of the Physics and Astronomy Division of the National 
Academy of Sciences of Ukraine. The authors acknowledge useful correspondence with 
Alexei Starobinsky. 



[1] D. Layzer, ApJ 138, 174L (1963); W. M. Irvine, PhD dissertation, Harvard University (1961). 
[2] N. A. Dmitriev and Ya. B. Zeldovich, Zh. Eksp. Teor. Fiz. 45, 1150 (1963) [Sov. Phys. JETP 
18, 793 (1963)]. 

[3] P. J. E. Peebles, The Large-Scale Structure of the Universe, Princeton University Press, Prince- 
ton (1980). 

[4] P. J. E. Peebles, Principles of Physical Cosmology, Princeton University Press, Princeton 
(1993). 

[5] F. Zwicky, Helv. Phys. Acta 6, 110 (1933). 

[6] L. Amendola and D. Tocchini-Valentini, Phys. Rev. D 66, 043528 (2002); L. Amendola and 

C. Quercellini, Phys. Rev. D 68, 023514 (2003). 
[7] M. Kesden and M. Kamionkowski, Phys. Rev. D74, 083007 (2006). 
[8] M. Fukugita and P. J. E. Peebles, ApJ 616, 643 (2004). 

[9] A. Nusser, S. S. Gubser and P. J. E. Peebles, Phys. Rev. D 71, 083505 (2005). 
[10] G. R. Farrar and P. J. E. Peebles, Astrophys. J. 604, 1 (2004). 

[11] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999); 83, 4690 (1999); P. Callin 

and F. Ravndal, Phys. Rev. D 70, 104009 (2004). 
[12] V. Sahni, A. Shafieloo and A. A. Starobinsky, Phys. Rev. D78, 103502 (2008). 
[13] V. Sahni and A. A. Starobinsky, IJMP(D) 15, 2105 (2006). 



